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Abstract 

The  rod-pinch  diode[  1 ,2]  is  a  cylindrical,  pinched-beam 
diode  being  developed  as  a  radiography  source[3].  The 
diode  consists  of  a  small  radius  anode  rod  extending 
through  the  hole  of  an  annular  cathode.  The  diode  has 
been  operated  at  1  to  5  MV  with  an  impedance  of  20  to  50 
Q,  a  FWHM  pulse  width  of  20  to  50  ns,  and  an  anode 
radius  as  small  as  0.25  cm[l-3].  The  diode  is  designed  to 
run  at  critical  current  so  that  electrons  emitted  from  the 
cathode  flow  axially  along  the  anode  rod  and  pinch 
radially  onto  the  rod  tip.  Typically,  ion  emission  from  the 
anode  is  required  for  propagation  of  the  pinch  along  the 
rod.  Without  ions,  the  pinch  would  occur  on  the  anode 
rod  just  downstream  of  the  cathode  disk.  In  order  to 
assure  that  a  given  diode  will  be  properly  designed  to  run 
at  critical  current  requires  a  detailed  knowledge  of  the 
diode  impedance  characteristics.  Initially,  a  laminar  flow 
model[4]  was  developed  to  describe  the  rod-pinch  diode. 
Although  this  model  provides  considerable  insight  into 
diode  operation,  PIC  simulations  show  that  the  electron 
flow  is  not  laminar[2].  The  model  of  [4]  was  extended  to 
include  transverse  electron  pressure  in  order  to  consider 
the  effects  of  nonlaminar  flow[5].  However,  a  form  for 
the  transverse  pressure  tensor  is  required  to  close  the 
equation  set  in  this  model  and  only  special  forms  of  the 
pressure  tensor  are  analytically  tractable.  Here,  a  Vlasov 
model  for  the  diode  electron  flow  is  developed  using  an 
electron  distribution  function  with  properties  that  are  well 
characterized  and  directly  related  to  a  rod-pinch  diode.  In 
this  model,  the  pressure  tensor  is  self-consistently  derived. 
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Figure  1.  Schematic  of  rod-pinch  diode. 

sheath  (rc  >  r  >  rs)  contains  ion  and  electrons,  while  the 
region  between  the  sheath  and  the  anode  (rs  >  r  >  ra) 
contains  only  ions.  At  z  =  l,  the  electrons  pinch  onto  the 
end  of  the  rod  (not  shown).  Because  t  »  rc,  d/dr  »  d/dz 
except  at  the  end  of  the  rod  where  the  pinch  occurs.  Thus, 
the  model  assumes  that  d/dz  =  0  and  canonical  axial 
momentum,  Pz  =  pz  -  eAz/c,  is  conserved,  where  pz  is  the 
axial  momentum,  Az(r)  is  the  vector  potential  with  Az(rc) 
=  0,  e  is  the  charge  on  an  electron  and  c  is  the  speed  of 
light.  Energy,  H  =  [m2c4  +  (pr2  +  pe2  +  pz2)c2]1/2  -  e<|>,  and 
canonical  angular  momentum,  Pe  =  rp@,  are  also 
conserved,  where  pr  and  pe  are  the  radial  and  azimuthal 
momenta,  4>(r)  is  the  electric  potential  with  4>(rc)  =  0,  and 
m  is  the  mass  of  the  electron.  For  electrons  born  on  the 
cathode  with  P0  =  0,  Pz  =  0,  and  H  =  me2,  the  electron 
distribution  function  is  written  as 


f  (p)  -  N 8( Pg  )8(P,  )8(H  -  me  2  )  ,  (1) 


I.  VLASOV  MODEL  FOR  THE 
ELECTRON  FLOW 

The  geometry  used  in  the  rod-pinch  diode  model  is 
illustrated  in  Fig.  1.  Ions  are  emitted  along  the  entire 
length  t  of  the  anode  rod  and  flow  radially  outward  while 
electrons  are  only  emitted  from  the  bottom  of  the  cathode 
at  radius  rc  and  flow  axially  in  a  sheath  above  the  anode. 
The  radius  of  the  anode  rod  is  ra  and  the  inner  edge  of  the 
anode  sheath  is  at  rs.  There  are  two  distinct  regions;  the 


where  N  is  a  constant.  Moments  of  the  distribution 
function  can  be  taken  to  derive  expressions  for  the 
electron  fluid  density  nc(r),  axial  velocity  Vez(r),  and 
pressure  Prr(r)  in  terms  of  Az(r)  and  4>(r).  In  evaluating 
these  integrals,  use  is  made  of  the  identity  5(g(x))  =  Z[5(x 
-  x;)/|dg(x)/dx|]  where  the  sum  is  over  the  zeroes  of  the 
function,  g(x;)  =  0.  These  expressions  are 

_N  ( e(j)  +  mc 2) 
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II.  EQUATIONS  FOR  (j>  AND  A- 


Prr=—[(e(j)  +  mc2)2 -e2A2z  -m2c4 ]1/2.  (4) 
cr 

As  a  self-consistency  check,  these  results  can  be  shown  to 
satisfy  the  momentum  transfer  equation.  All  other 
components  of  the  fluid  velocity  and  pressure  tensor  are 
identically  zero.  Note  that  the  pressure  Prr  given  in  Eq. 
(4)  vanishes  both  at  the  cathode  where  the  electron  have 
zero  kinetic  energy  (and  <|>  =  Az  =  0)  and  at  the  sheath 
edge  where  electrons  move  only  axially.  This  sheath 
condition,  Prr(rs)  =  0,  will  be  needed  to  obtain  a  solution. 
As  expected,  Pn(r)  is  found  to  have  a  maximum  between 
the  cathode  and  the  sheath  edge.  Also  note  that  it  follows 
from  the  properties  of  Prr  that  nc  is  singular  at  both  the 
cathode  and  sheath  edge  and  has  a  minimum  between  the 
cathode  and  sheath  edge.  Expansion  techniques  to  obtain 
approximate  analytic  solutions  at  these  singularities  are 
used  to  aid  in  obtaining  the  numerical  solution.  Finally, 
Vez  vanishes  at  the  cathode  where  Az  =  0  and,  as  derived 
from  the  sheath  condition,  increases  to  (ys2  -  l)1/2c/ys  at  the 
sheath  edge  where  ys  =  1  +  e(|)(rs)/mc2  is  the  electron 
relativistic  factor  at  the  sheath  edge. 


1  d  {ndy\-  JecY 

RdRy  a RJ  R(y2  -  S2  - 1  )I/2 


(5) 


1  d  nds\  JecS 

- R —  = - ; - - ttt  .  (6) 

RdR{  dRJ  R(y2  -  S2  - 1  )1/2 

where  R  =  r/rc,  y(R)  =  1  +  e(|)(R)/mc2,  S(R)  =  eAz(R)/mc2, 
Jec  =  47ie2Nrc/me'  is  the  normalized  axial  electron  current 
density  at  rc,  Jic  =  (47tejicrc2/me3)(M/2mZ)1/2  is  the 
normalized  radial  ion  current  density  at  rc,  jic  is  the  radial 
ion  current  density  at  rc,  and  M  and  Z,  are  the  ion  mass 
and  charge  state.  The  ion  contribution  on  the  right  hand 
side  of  Eq.  (5)  is  the  same  as  derived  previously[4,5].  In 
the  region  below  the  sheath  (i.e.,  rs  >  r  >  ra),  the  same 
equations  apply  but  with  Jec  =  0.  Since  d/dz  =  0  was 
assumed,  the  local  value  of  the  ion  current  (through  the 
boundary  conditions)  determines  the  axial  location  z 
where  the  solution  applies.  For  example,  the  full  ion 


Figure  2.  Solution  for  case  with  ya  =  4.9139  (i.e.,  V  =  2  MV),  Ra  =  ra/rc  =  0.1,  and  l/rc  =  2  where  Rs  =  rs/rc  is  found  to  be 
0.155.  Flere,  Ne  =  47irc2eenc/IA,  Ni  =  47irc2eZenc/IA,  and  P  =  47iercPrr/meIA. 
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current  flows  in  the  rod  at  the  cathode  location  (z  =  0)  and 
no  ion  current  flows  in  the  rod  at  the  end  of  the  rod  (z  = 
£).  A  series  of  such  local  solutions  provides  a  global 
solution.  Here,  only  solutions  at  the  cathode  location  are 
considered,  so  that  z  =  0  and  I;  =  2jircfjic,  where  it  is 
assumed  for  simplicity  that  j;c  does  not  vary  with  z. 
Solutions  are  obtained  by  integrating  Eqs.  (5)  and  (6) 
from  the  outer  boundary  in  each  region  (i.e.,  from  rc  and 
ra)  to  the  sheath  where  R  =  rs/rc  =  Rs.  Three  boundary 
conditions  are  known  at  the  cathode  where  R  =  1  and  two 
are  known  at  the  anode  where  R  =  ra/rc  =  Ra.  These  are 
y(l)  =  1,  y'(l)  =  0,  S(l)  =  0,  y(Ra)  =  1  +  eV/mc2,  and 


y'(Ra)  =  0  where  V  =  c|j(Ra).  There  are  five  unknowns  Jec, 
Rs,  S'(l),  S(Ra),  and  S'(Ra)  which  must  be  guessed  to 
carry  out  the  integrations.  The  three  unknown  boundary 
conditions  can  be  expressed  as  -S'(l)  =  eB(rc)rc/mc2  = 
2Ii/IA,  S(Ra)  =  eAz(ra)/mc2,  and  -S'(Ra)  =  eB(ra)rc/mc2  = 
2(1;  +  Ie)/IARa.  Here,  I;  and  Ie  are  the  ion  and  electron 
currents  and  IA  =  mc3/e  =17  kA  is  related  to  the  Alfven 
current.  Note  that  Jic,  which  is  needed  to  integrate  Eq.  (5), 
can  be  written  in  terms  of  S'(l)  through  J;c  = 
(2I;rc/fIA)(M/2Z;m)1/2.  Valid  values  for  these  five 
unknowns  are  determined  by  five  conditions  that  need  to 
be  satisfied.  These  include  matching  y(Rs),  y'(Rs),  S(Rs) 


Figure  3.  Comparisons  of  numerical  solutions  with  PIC  simulation  results,  other  theoretical  predictions,  and  Sabre  data. 
In  (a)  and  (d),  V  =  2  MV  and  Urc  =  2;  in  (b)  and  (e),  V  =  2  MV  and  rc/ra  =  10;  and  in  (c)  and  (f),  f/rc  =  2  and  rc/ra  =  10. 
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and  S'(Rs)  from  the  two  solutions  at  the  sheath  and 
satisfying  the  sheath  condition  [y2(Rs)  -  S2(Rs)  -  1]  =  0 
(from  Pn-(RS)  =  0).  Thus,  the  problem  is  well  posed  and 
the  solution  provides  values  for  the  five  unknowns  or 
equivalently  for  Jec,  Rs,  Az(ra),  I,,  and  Ie. 


III.  RESULTS  AND  DISCUSSION 

Results  from  the  numerical  solution  of  Eqs.  (5)  and  (6) 
are  shown  in  Fig.  2  for  the  case  where  ya  =  4.9139,  rc/ra  = 
10,  and  H rc  =  2.  Here,  the  normalized  densities  are  Nc  = 
47irc2ecne/IA  and  N  =  47rrc2eZcn;/IA,  and  the  normalized 
pressure  is  P  =  47irl2ePrr/mcIA.  Note  that,  as  predicted,  the 
electron  density  diverges  at  the  cathode  and  at  the  edge  of 
the  sheath  and  the  electron  pressure  goes  to  zero  at  the 
cathode  and  the  edge  of  the  sheath.  Also,  in  the  region 
between  the  anode  and  the  sheath  where  there  are  only 
ions  the  electric  field  is  highly  peaked.  Results  obtained 
under  various  conditions  are  compared  with  PIC 
simulation  results[2],  other  theoretical  predictions[2],  and 
data  from  experiments  on  the  Sabre  accelerator  at  Sandia 
National  Laboratories[3,6]  in  Fig.  3.  Plots  of  a,  Z,  f,  and 
Ie  are  shown  over  a  range  of  values  for  rc/ra,  f/rc,  and  ya. 
The  electrical  performance  of  the  diode  is  characterized 
by  its  impedance  Z  =  V/I,  where  I  =  Ie  +  I;  is  the  total 
diode  current.  Values  of  Z  obtained  from  the  numerical 
solutions  compare  reasonably  well  with  the  data  obtained 
on  Sabre  as  shown  in  Figs.  3a,  3b,  and  3c.  The  diode 
current  has  been  predicted[2]  to  vary  as  I  =  alcrjt  = 
8500a(y2  -  l)1/2/ln(rc/ra)  where  a  is  a  scaling  factor. 
Values  of  a  =  (Ie  +  Ii)/Icrit  from  the  numerical  solutions  are 
close  in  value  but  scale  somewhat  differently  from  the 
factor  aPIC  found  in  PIC  simulations  as  illustrated  in  Figs. 
3a,  3b,  and  3c.  Similarly,  the  ion  current  has  been 
theoretically  predicted[2]  to  scale  roughly  as  I;th  = 
1.6Ic[Zjm(y  +  l)/2M]1/2f/(rc  -  ra)  for  t »  (rc  -  ra).  Without 
enhancement  from  the  pinched  electron  flow,  the  space- 
charge-limited  ion  current  IiSCL  =  7.33  x  10"3(ya  - 
1  )3  2Cp2rc  would  flow.  Here,  (32  is  given  in  [7]  and  the 
ions  are  protons.  Values  of  Ie  and  I;  from  the  numerical 
solutions  are  compared  with  l/h  (using  Ie  from  the 
numerical  solution  in  the  above  equation  for  I;th)  and  I;SCL 
in  Figs.  3d,  3e,  and  3f.  As  expected,  values  of  I;  found 
from  the  numerical  solutions  are  much  larger  than  I;SCL. 
The  agreement  with  I*  is  reasonable  in  absolute  value  but 
the  scaling  is  somewhat  different. 


IV.  SUMMARY 

A  Vlasov  model  for  the  diode  electron  flow  in  a  rod- 
pinch  diode  is  developed  that  is  derived  from  a 
distribution  function  for  the  electrons  with  properties  that 
are  well  characterized  and  directly  related  to  the  rod-pinch 
diode.  Non-laminar  flow  is  self-consistently  included  and 
results  in  a  transverse  radial  pressure  which  is  defined  in 


Eq.  (4).  Results  from  the  model  agree  reasonable  well  in 
magnitude  with  the  data  and  results  from  PIC  simulations 
but  differences  in  scaling  are  observed.  In  the  model,  all 
electrons  are  born  with  only  their  rest  energy  and  zero 
canonical  angular  momentum  as  would  be  the  case  in  an 
experiment  where  all  electrons  are  born  at  rest  on  the 
cathode  surface.  However,  electrons  in  the  model  are 
restricted  to  have  Pz  =  0,  which  is  only  true  for  electrons 
bom  on  the  bottom  of  the  cathode.  Electrons  born  on  the 
upstream  and  downstream  faces  of  the  cathode  would 
have  different  Pz.  It  is  also  assumed  that  d/dz  =  0,  so  that 
Pz  is  conserved.  This  assumption  is  also  not  strictly 
correct.  It  is  possible  that  the  true  2D  nature  of  the 
problem  including  emission  from  the  sides  of  the  cathode 
(i.e.,  Pz  0  for  all  electrons)  and  axial  gradients  (i.e.,  d/dz 
^  0)  in  the  fields  near  the  cathode  and  the  pinch  at  the  end 
of  the  rod  can  explain  the  observed  differences  in  scaling. 
The  simulation  results  displayed  in  Fig.  3  are  from  cases 
where  the  simulation  geometry  is  set  up  to  model  the 
actual  experimental  geometry.  Simulations  will  be  run  in 
geometries  that  more  closely  resemble  the  case  described 
by  the  model  to  study  these  issues. 
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